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1 .  INTRODUCTION 


This  is  the  Final  Technical  Report  on  a  program  conducted 
by  Nielsen  Engineering  &  Research,  Inc.  (NEAR)  for  the  Office  of 
Naval  Research  to  investigate  various  aspects  of  the  structure 
and  computation  of  incompressible  turbulence.  Studies  have  been 
conducted  over  a  five-year  period  in  three  main  areas:  (1) 
testing  of  subgrid-scale  (SGS)  models  for  use  in  large-eddy 
simulation  (LES)  of  homogeneous  turbulent  flows;  (2)  investiga¬ 
tions  into  the  structure  of  homogeneous  shear  flow;  and  (3) 
developing  a  model  of  turbulence  near  a  wall  from  solutions  of 
the  Navier-Stokes  equations . 

In  t^he  first  two  of  these  studies,  we  have  been  among  the 
first  investigators  to  use  a  technique  which  will  assume  increas¬ 
ing  importance  in  the  years  to  come:  the  use  of  full  numerical 
solutions  of  the  Navier-Stokes  equations  as  "data"  for  model 
testing,  determination  of  model  parameters,  and  studies  of  the 
basic  flow  dynamical  structure.  For  simple  enough  flows  and  if 
done  with  sufficient  care,  these  full  simulations  (which  use  no 
models)  can  be  considered  exact  representations.  As  such,  they 
provide  information  on  quantities  of  interest  in  unprecedented 
detail,  and  allow  the  study  of  quantities  that  currently  cannot 
be  measured  experimentally.  In  our  studies,  we  have  used  some 
of  the  full  simulations  of  homogeneous  flows  done  by  Dr.  Rogallo 
of  the  NASA/Ames  Research  Center  (ref.  1)  . 

The  first  study  has  involved  the  use  of  certain  of  these 
simulations  as  "data"  to  evaluate  the  models  used  in  LES.  In 
LES,  the  largest  structures  in  the  flow  are  computed  directly 
using  time-dependent  three-dimensional  equations  derived  from 
the  Navier-Stokes  equations  by  the  application  of  a  spatial 
filter.  The  filtering  introduces  terms  involving  the  small-scale 
motions  which  must  be  modeled  (the  LES  version  of  the  turbulence 
closure  problem).  In  our  work,  we  have  tested  various  SGS  models 
by  treating  the  results  of  Rogallo 's  full  simulations  as  an 
experimentalist  would  treat  measurements.  The  "data"  from  a 


given  time  step  in  the  full  simulation  are  filtered  to  give  the 
large-scale  component  of  the  flow  field,  and  the  small-scale 
field  is  obtained  by  subtraction.  The  terms  which  need  to  be 
modeled  are  then  computed,  as  well  as  the  model  prediction  of 
these  terms  using  the  large-scale  field,  and  a  direct  comparison 
between  the  exact  and  modeled  values  is  then  made.  The  use  of 
this  approach  to  evaluate  various  forms  of  SGS  models  is  iso¬ 
tropic  turbulence  and  in  turbulence  interacting  with  mild  strain 
was  reported  in  references  2-4.  The  extension  of  this  work  to 
cases  involving  strong  strain  or  shear  is  contained  in  refer¬ 
ence  5  and  in  Section  2  of  this  report. 

The  second  study  undertaken  under  this  research  program 
also  uses  the  full-simulation  results  of  Rogallo  as  "data."  In 
this  investigation,  attempts  have  been  made  to  determine  the 
mechanisms  underlying  the  production  of  turbulence  energy  and 
growth  of  length  scales  in  homogeneous  shear  flow  using  computer¬ 
generated  displays  of  selected  quantities  as  a  form  of  "flow 
visualization."  Progress  achieved  to  date  in  this  study  is  dis¬ 
cussed  in  Section  3.  Because  most  of  the  work  to  this  point  in 
this  study  has  involved  developing  the  computer-graphic  tools 
required,  only  preliminary  results  are  available  at  this  time; 
this  work  is  scheduled  to  be  continued  in  another  program  under 
ONR  sponsorship. 

The  aim  of  the  third  study  is  the  development  of  a 
computational  model  of  the  viscous  sublayer  of  wall-bounded 
incompressible  turbulent  flow.  In  this  work,  the  time-dependent 
Navier-Stokes  equations  (simplified  in  accord  with  the  experi¬ 
mental  observations  of  the  sublayer)  are  used  to  compute  the 
characteristics  in  the  sublayer  subject  to  space-  and  time- 
dependent  boundary  conditions  specified  at  the  outer  edge  of  the 
sublayer.  The  essence  of  this  work  is  the  formulation  of  these 
boundary  conditions;  experimental  measurements  have  been  used 
to  guide  their  construction,  and  considerable  success  has  been 
achieved  with  regard  to  the  agreement  with  data  of  the  resulting 
distributions  throughout  the  sublayer  of  most  of  the  important 
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turbulence  quantities.  This  study  is  reported  in  reference  6 
and  Section  4  and  will  also  be  continued  under  ONR  sponsorship. 


2.  TESTS  OF  SUBGRID-SCALE  MODELS  IN 
STRAINED  TURBULENCE 


In  references  1  and  5,  it  is  shown  that  Rogallo's  full 
simulations  of  strained  and  sheared  turbulence  are  representa¬ 
tive  of  physical  flows,  as  long  as  the  simulations  are  not 
carried  too  far  in  time.  The  reasons  for  this  are  related  to 
the  fact  that  the  computations  are  done  using  transformations 
of  the  independent  and  dependent  variables,  so  that  the  physical 
coordinate  system  distorts  as  the  computation  proceeds.  In 
Rogallo's  simulation,  the  flow  is  divided  into  an  imposed  mean 
flow  (U^)  and  the  turbulence  (u^) .  The  mean  flow  is  prescribed 
to  be  of  the  form: 


U. 

l 


r . .  x, 
ik  k 


(1) 


and  the  solution  of  the  equations  for  u.  (derived  from  the 
Navier-Stokes  equations)  is  done  after  transformation  of  the 
independent  and  dependent  variables : 


£i  =  Bikxk  Khere  Eik  +  Eijrjk  =  0 


v^  =  (no  summation) 


(2) 

(3) 


Most  of  the  solutions  were  done  using  64  mesh  points  in  each 
direction,  although  a  few  calculations  (discussed  later)  were 
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done  using  a  128  mesh;  m  either  case,  the  calculations  were 
done  entirely  in  Fourier  space. 

For  pure  steady  plane  strain  (T^  =  T,  1^2  ~  •  aH  other 

I\j  =  0)  the  result  of  the  transformation  (2)  is  that  the  compu¬ 
tational  region  ultimately  becomes  very  short  in  the  2-direction 

and  the  calculation  must  be  stopped  when  the  total  strain  ratio 
Ft 

e  reaches  a  value  near  4.  Continuing  the  calculations  beyond 
this  point  results  in  the  calculated  flow  becoming  correlated 
in  the  compressed  direction,  which  is  of  course  unphysical. 


Additionally,  the  distortion  produces  a  grid  for  which  all  the 
wave  numbers  in  the  stretched  direction  represent  relatively  large 
scales,  the  dissipation  range  is  not  properly  captured,  and  proper 
energy  transfer  is  impeded.  While  strain  ratios  greater  than  four 
are  often  achieved  experimentally  (see,  for  example,  ref.  7),  a 
value  of  four  is  large  enough  to  be  useful.  In  the  case  of  steady 
shear  (T.^  =  °ther  =  0),  the  region  shears  and  must 

periodically  be  readjusted.  The  readjustment  causes  some  loss 
of  information  and  requires  that  the  results  be  treated  carefully. 

In  most  of  the  calculations  used  as  a  basis  for  model  test¬ 
ing,  plane  strain  was  the  force  driving  the  turbulence  away  from 
isotropy.  A  number  of  shear  cases  were  also  run  and  no  signifi¬ 
cant  qualitative  difference  between  shear  and  strain  was  found. 
Homogeneous  shear  of  strength  T  can  be  regarded  as  a  combination 
of  plane  strain  of  strength  T/2  at  45°  and  rotation  at  angular 
velocity  T/2.  Since  the  effect  of  rotation  is  to  reduce  the 
production  rate  in  the  turbulence  (cf .  ref.  8) ,  we  expect  the 
shear  of  strength  T  will  produce  effects  similar  to,  but  slightly 
weaker  than  strain  of  strength  T/2;  this  is  borne  out  by  the 
calculations . 

Since  the  earlier  results  of  references  3  and  4  indicated 

that  flows  with  low  strain  rates  are  not  significantly  different 

from  isotropic  turbulence  and  because  most  technological  flows 

experience  much  higher  strain  rates,  we  decided  to  proceed  by 

successively  doubling  the  strain  rate.  In  the  end,  the  basic 

matrix  contained  four  strain  rates:  1.33,  2.65,  5.30,  and  10.61 

sec  For  each  strain  rate,  the  calculation  was  carried  to 

Ft 

total  strain  e  of  1.3,  2,  and  4.  The  initial  Reynolds  number 
was  kept  at  40  in  all  cases  (the  maximum  available  using 
Rogallo’s  simulation).  Thus  a  total  of  twelve  cases  were  avail¬ 
able  for  analysis;  of  these,  ten  were  actually  analyzed.  Several 
additional  runs  were  made  for  special  purposes.  It  is  no  sur¬ 
prise  that  the  most  interesting  cases  turned  out  to  be  the  ones 
with  the  highest  values  of  both  strain  rate  and  total  strain. 


To  use  the  results  of  the  full  simulations  in  our  studies 
of  subgrid-scale  modeling,  the  following  definitions  are  required, 
The  large-scale  (filtered)  velocity  field  is  obtained  from 


=i-J 


G (r  -  r ' ) u± (r ' ) dr ' 


where  G  is  a  filter  function  which  is  usually  chosen  to  be  a 

2  -3/2  22 

Gaussian,  (ttA  /6)  exp(-6r  /A  )  .  In  this  work  we  use  this 

3  3 

filter  with  A  /A  =  2;  A  is  the  spacing  of  the  grid  on  which 

3  C  C 

the  modeling  is  done.  The  subgrid-scale  velocity  field  is  then 
defined  by  uj  =  u^  -  iL  and  the  subgrid  Reynolds  stress  is 
defined  by 

Tij =  1 5ij(2“l^k +  “Wi*  •  “I5! 


-  u.u!  -  u ! u !  =  =•  6  .  .  R,  ,  -  R.  . 

ID  ID  3  id  kk  id 

We  will  also  look  separately  at  the  term 

t .  .  =  4  ^ •  u ' u '  -  uTuT  =  4  6 •  •  rLi,  -  R-  •  (6) 

ID  3  id  k  k  id  3  id  Kk  id 

which  has  been  called  the  backscatter  term  (ref.  9) . 

One  can  derive  equations  for  the  subgrid-scale  Reynolds 
stresses;  these  are  similar  to  the  equations  for  the  usual 
time-averaged  Reynolds  stresses  that  are  commonly  used  in  turbu¬ 
lence  modeling.  From  these  equations  one  can  derive  an  equation 
for  the  subgrid-scale  kinetic  energy  <iLu|  +  u^u^/2>.  The  pro¬ 
duction  term  in  this  equation  is  <TT~W7/9Xj>  and  is  the  energy 
transfer  from  the  large  to  small  scales.  We  have  computed  this 
term  and  additionally  have  separated  the  contributions  from 

and  t. ..  These  results  are  shown  in  table  1  for  two  flow  fields 

ID  _ = - 

with  a  strain  ratio  of  4.  Examination  of  <T^9u^/9Xj>  reveals 

that  for  the  smaller  of  the  two  strain  rates  shown,  the  energy 

transfer  is  in  the  normal  direction  (to  the  small  scales) ,  but 

this  reverses  at  the  higher  strain  rate.  Comparison  with  energy 

transfer  via  <f . . 9u . /9x . >  shows  that  this  reversal  is  due  to  the 
ID  i  D 

backscatter  term,  which  for  these  two  flows  is  transferring 
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energy  from  the  small  scales  to  the  large  ones.  The  energy 
transfer  due  to  the  "cross"  or  "outscatter"  terms  (the  differ¬ 
ence  between  <t. .3u./3x->  and  <f. .3u./3x->)  is  seen  to  be  in  the 
13  1  j  13  i  j 

normal  direction  for  both  flows  of  table  1  but  is  overwhelmed  by 
the  backscatter  for  r  =  10.61  sec  ^ ,  resulting  in  a  net  transfer 
of  energy  to  the  large  scales.  The  behavior  of  the  backscatter 
term  is  shown  in  figure  1  over  the  whole  range  of  strain  ratios 
covered  for  the  two  strain  rates  of  table  1.  Energy  transfer 
from  the  small  scales  to  the  large  scales  seems  to  be  associated 
with  high  strain  rates  which  act  on  the  flow  for  a  long  time. 

The  existence  of  energy  transfers  from  the  subgrid-scale  to 
the  large  scales  has  a  considerable  effect  on  the  validity  of 
the  subgrid-scale  models  as  will  be  shown  in  the  next  section. 
Repeat  calculations  with  different  initial  conditions  (but  the 
same  initial  spectrum)  have  shown  that  this  effect  is  not  the 
consequence  of  peculiar  initial  conditions.  To  eliminate  the 
possibility  that  this  phenomenon  is  a  computational  effect,  the 
highest  strain-rate  case  was  recomputed  more  accurately  by 
Rogallo  using  a  1283  mesh.  The  results  for  this  case  are  also 
shown  in  table  1:  it  is  clear  that  the  energy  transfer  to  the 
large  scales  persists  in  this  more  accurate  calculation,  and  is 
in  fact  more  vigorous.  Limited  tests  were  also  made  with  an 
alternate  filter  definition  (a  sharp  cut-off  filter  in  wave 
space) .  These  suggest  that  although  the  amount  of  energy  trans¬ 
fer  between  scales  is  obviously  dependent  on  the  definition  of 
the  filter  used  and  on  its  width,  the  gross  direction  of  the 
energy  transfer  is  not.  In  pseudo-spectral  calculations  of  homo¬ 
geneous  shear  flows,  Shirani,  et  al.  (ref.  10)  have  observed  a 
similar  phenomenon;  in  their  cases,  there  was  energy  transfer 
from  both  small  and  large  scales  to  the  medium-sized  eddies. 

These  results,  which  are  not  in  accord  with  the  usual  concept  of 
the  energy  "cascade"  towards  ever-smaller  scales,  may  be  of  con¬ 
siderable  importance  in  turbulence  modeling. 

The  methodology  used  to  test  the  models,  which  is  an  exten¬ 
sion  of  that  of  reference  11,  is  now  briefly  described.  From  the 
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velocity  fields  described  previously,  which  we  accept  as  exact 

realizations  of  a  turbulent  flow  field,  we  can  compute  (or 

t. .)  on  a  coarse  grid  (A  =  4A)  overlaid  on  the  computational 
1 J  c 

grid.  Furthermore,  we  can  compute  model  versions  of  these  quan¬ 
tities  which  are  based  on  parameters  of  the  large-scale  field. 

A  direct  test  of  the  model  is  made  by  computing  the  correlation 
coefficient  n  on  the  coarse  grid  between  the  exact  subgrid-scale 
Reynolds  stress  and  the  model  for  the  stress  itself  (the  "tensor 
level") ,  its  divergence  (the  "vector  level") ,  and  its  associated 
energy  transfer  (the  "scalar  level") .  For  models  which  depend 
linearly  on  a  single  constant,  the  correlation  coefficients  are 
independent  of  the  model  constant  and  thus  provide  an  unbiased 
test  of  the  model.  The  constant  can  then  be  found  by  least- 
squares  fitting  but  it  was  found  that  matching  the  mean-square 
values  of  the  model  and  the  exact  result  produces  more  satisfac¬ 
tory  results  in  most  cases. 


2.1  Eddy  Viscosity  Models 


Almost  all  large-eddy  simulation  work  to  date  has  used  eddy 
viscosity  models  in  which  the  effect  of  the  subgrid-scale  on  the 
large-scale  turbulence  is  represented  as: 


T 


ij 


(7) 


where  vT  is  the  eddy  viscosity.  A  number  of  formulations  of  the 
eddy  viscosity  have  been  given  by  various  authors  but  experience 
with  them  and  the  direct  tests  of  references  4,  5,  and  11  have 
shown  that  there  is  in  fact  little  to  choose  among  the  various 
possibilities.  Although  eddy  viscosity  models  are  superior  in 
some  ways  to  other  models  that  have  been  suggested,  they  are  far 
from  ideal.  We  will  need  to  make  do  with  them  until  more  accur¬ 
ate  models  are  proposed  and  tested.  Thus  the  first  task  is  to 
extend  the  method  of  references  4  and  5  to  the  strain  and  shear 
flows  described  above.  Discussion  of  other  kinds  of  models  will 
follow. 
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We  will  focus  on  a  particular  formulation  of  the  eddy  vis¬ 
cosity,  the  Smagorinsky  model: 

VT  =  (CsAa)2^l  (8) 

where  |s|  is  the  mean  strain  of  the  large-scale  field  and  is 
the  filter  width. 

In  our  previous  work,  for  isotropic  decay  the  correlation 

coefficient  and  model  parameters  were  found  to  be  functions  of 

a  single  nondimensional  parameter  —  the  subgrid-scale  Reynolds 

number  R  ;  time  was  found  not  to  be  important  in  this  flow, 
sgs 

The  model  parameter,  but  not  the  correlation  coefficient,  was 
also  found  to  depend  somewhat  on  the  numerical  method  used  for 
differentiation  in  the  model  calculations  but  this  is  a  compu¬ 
tational  effect  rather  than  a  physical  one.  The  presence  of 
strain  means  that  the  strain  rate,  F,  is  also  an  important  inde¬ 
pendent  parameter.  Since  it  is  dimensionally  inverse  time,  it 
can  be  nondimensionalized  using  a  time  scale  and  there  are  two 
obvious  choices.  The  first  is  the  actual  time  and  the  nondimen- 

rt 

sional  quantity  formed  using  it  is  Tt;  e  ,  the  total  strain,  is 
a  parameter  of  obvious  importance  in  these  flows.  The  second 
time  scale  is  one  characteristic  of  the  flow  itself.  The  best 
choice  is  td,  the  dissipation  time  scale,  the  ratio  of  the  turbu¬ 
lence  kinetic  energy  to  the  dissipation  rate.  Since  the  Reynolds 
number  is  fixed  in  this  work,  we  have  two  parameters  with  which 
to  correlate  the  results. 

In  the  Smagorinsky  model  it  is  not  clear  whether  the  imposed 

mean  strain  rate  should  be  included  in  the  model  (in  vT  and  S^) 

or  not.  On  the  one  hand,  one  can  argue  that  in  a  large-eddy 

simulation  it  is  difficult  to  distinguish  the  mean  strain  from 

the  turbulent  component.  On  the  other  hand,  if  we  consider  the 

homogeneous  strain  flows  with  very  large  imposed  strain,  rapid 

distortion  theory  is  approximately  correct  and  the  state  of  the 

Ft 

turbulence  will  depend  primarily  on  the  total  strain  e  .  The 
model  parameter  will  then  be  inversely  proportional  to  the 


imposed  strain  at  large  strain  rates.  We  have  therefore  decided 
to  test  the  model  both  with  and  without  the  mean  strain  included. 


The  scalar  level  correlation  coefficients*  for  the 

Smagorinsky  model  are  plotted  as  functions  of  the  total  strain 

ratio  in  figure  2.  The  results  with  and  without  mean  strain  in 

the  model  are  essentially  identical  so  that  there  is  nothing  to 

choose  between  them  on  this  basis.  The  value  for  isotropic  decay 

is  shown  for  comparison.  We  see  that  the  validity  of  the  model 

becomes  dubious  at  large  total  strains.  The  negative  value  of 

the  correlation  coefficient  at  large  strains  reflects  in  part 

the  reversed  direction  of  energy  transfer  remarked  on  earlier. 

These  results  were  obtained  using  t . .  rather  than  t . .  but  calcu- 

il  ID 

lations  obtained  using  show  no  substantial  differences. 

The  model  parameter  for  the  Smagorinsky  model  with  the  mean 
strain  included  is  shown  in  figure  3;  figure  4  shows  the  same 
quantity  without  the  mean  strain.  Again  the  parameter  calculated 
for  isotropic  decay  is  shown  for  comparison.  The  reduction  in 
the  model  parameter  at  large  strain  rates  when  the  mean  strain  is 
included  in  the  model  is  the  effect  anticipated  above.  When  the 
mean  strain  is  excluded  the  model  parameter  is  very  nearly  inde¬ 
pendent  of  the  imposed  mean  strain  rate.  We  thus  conclude  that 
it  is  better  not  to  include  the  mean  strain  in  the  model  but,  as 
remarked  earlier,  this  may  be  very  difficult  to  do  in  a  large- 
eddy  simulation  of  an  inhomogeneous  shear  flow.  Perhaps  a  model 
similar  to  the  two-term  model  proposed  in  reference  12  would  have 
the  kind  of  properties  that  we  desire. 

The  behavior  of  the  correlation  coefficients  is  extremely 
disappointing  and  means  that  eddy  viscosity  models  cannot  be 


♦These  correlation  coefficients  were  calculated  without  re¬ 
moving  the  mean  of  the  exact  or  modeled  scalar;  it  is  better 
statistical  practice  to  correlate  the  exact  and  modeled  values 
of  Ui3fij/3xj-<ui3fij/3xj>.  Doing  this  results  in  correlation 
coefficients  which  are  lower  than  those  shown  by  as  much  as  0.16, 
but  the  conclusions  demonstrated  by  figure  2  are  unaltered. 


trusted  in  flows  with  appreciable  mean  strain.  (Although  we 
have  studied  only  the  Smagorinsky  model  here,  similar  results 
have  been  obtained  for  the  other  eddy  viscosity  models.) 

Another  way  to  see  the  inability  of  the  model  to  represent 
the  subgrid-scale  stresses  is  shown  in  figures  5  and  6.  These 
"scatter  plots"  show  the  exact  vs.  the  modeled  values  of 
ui3Tij/3Xj  on  the  coarse  grid  and  clearly  indicate  that  the 
problem  is  truly  the  lack  of  correlation  between  the  model  and 
the  exact  result;  no  simple  change  in  the  model  will  help.  On 
this  type  of  plot,  exact  agreement  between  the  modeled  and  exact 
quantities  would  be  represented  by  the  plotted  symbols  falling 
on  a  line  bisecting  the  first  and  third  quadrants.  This  is 
obviously  far  from  the  case  here.  (These  plots  were  produced 
on  microfiche  by  the  computer.  In  the  shaded  area,  the  symbols 
fell  so  close  together  that  resolution  of  individual  points  was 
lost.)  The  model  obviously  does  not  represent  the  subgrid-scale 
turbulence  accurately  even  at  low  strain  rates;  the  correlation 
coefficient  of  .45  is  actually  not  very  high.  We  also  made  plots 
showing  the  relative  directions  of  the  principal  axes  of  the 
model  and  the  exact  stress  and  found,  as  anticipated,  that  the 
two  sets  of  axes  appear  to  be  essentially  random  with  respect  to 
each  other. 

These  results  may  explain  why  the  Smagorinsky  model  has  not 
been  very  effective  in  computing  inhomogeneous  shear  flows, 
especially  the  channel  flow  which  has  been  computed  by  a  number 
of  authors  (refs.  12-15) .  There  is  clearly  a  need  for  improve¬ 
ment  in  the  models  used.  Recently,  a  new  concept  in  modeling 
was  introduced  (ref.  16) ,  and  we  have  performed  some  preliminary 
tests  on  this  model.  Some  of  these  results  are  presented  in  the 
next  section;  more  comprehensive  evaluation  is  to  be  reported 
in  reference  17. 


2.2  New  Models 

In  reference  16,  the  idea  is  advanced  that  the  major  inter¬ 
action  between  the  large-  and  small-scale  components  is  between 
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the  smallest  scales  of  the  large-scale  field  and  the  largest 
scales  of  the  small-scale  field  ("scale  similarity").  Since 
represents  the  large-scale  component  of  the  field,  applying  the 
filter  again  produced  a  field  (u^)  still  richer  in  the  largest 
scales,  and 

ui  =  -  ui  (9) 


is  a  field  which  contains  the  smallest  scales  of  the  large-scale 
component.  A  simple  version  of  a  model  which  incorporates  the 
notion  of  scale  similarity  is  then 


: .  .  =  C  ( =-6  .  .  u,  u,  -  u .  u  . ) 
i]  3  i]  k  k  i  j 


(10) 


Some  exploratory  results  using  equation  (10)  are  shown  in  tables  2 

and  3.  In  table  2,  the  correlation  coefficients  for  this  model 

in  homogeneous  isotropic  turbulence  at  the  tensor  (comparison  with 

3t  .  . 

f .  .)  ,  vector  (comparison  with  1 3-)  ,  and  scalar  (comparison  with 

1 J  ox  • 

- 

ui  9x  )  levels  are  given.  Values  for  the  Smagorinsky  model  are 

given  for  comparison  and  it  is  clear  that  the  model  of  equa¬ 
tion  (10)  has  considerable  promise.  An  even  more  impressive 
comparison  is  given  in  table  3  which  presents  the  analogous 
quantities  for  strained  turbulence.  In  these  tables  the  corre¬ 
lation  coefficients  were  calculated  after  removing  the  means  (in 
contrast  to  those  of  fig.  2) .  A  typical  scatter  plot  is  given  in 
figure  7;  the  improvement  over  the  previous  results  is  obvious. 


A  version  of  the  scale-similarity  model  which  is  seemingly 
even  more  accurate  is  described  in  reference  16: 


a-3 


•  c[j(ukuk  -  'k“k)Sij  -  (uiuj  -  ”iV]  (11) 


In  that  reference,  correlation  coefficients  for  isotropic  and 
sheared  turbulence  are  shown  which  are  even  higher  than  those 
for  the  model  of  equation  (10)  in  tables  2  and  3. 


The  apparent  success  of  these  models  led  to  a  test  in  a 
large-eddy  simulation  by  Bardina,  et  al.  and  they  found  that 


these  models  are  not  dissipative.  As  a  result,  a  search  was  made 
for  models  which  contained  the  best  features  of  the  new  models 
and  the  Smagorinsky  model.  Among  the  ideas  that  we  tried  was  a 
Smagorinsky  model  with  the  large-scale  velocity  u^  replaced  by 
u^;  the  results  were  found  to  be  very  similar  to  those  for  the 
Smagorinsky  model.  Another  idea  was  to  model  directly 

in  terms  of  CL  ;  this  model  was  found  to  have  almost  no  validity 
at  all. 

The  suggestion  was  then  made  that  a  model  which  is  a  linear 
combination  of  the  Smagorinsky  model  and  equation  (11)  should 
have  the  best  features  of  both  models.  Some  promising  preliminary 
tests  of  this  concept  have  been  made  and  investigation  of  com¬ 
bined  models  of  this  type  has  been  continued  (ref.  17) . 

3.  STUDIES  OF  STRUCTURE  OF  HOMOGENEOUS 
SHEAR  FLOWS 

In  the  past  decade,  vortex  pairing  has  been  put  forward  as 
a  primary  mechanism  for  the  growth  of  free-shear  flows  (particu¬ 
larly  in  the  mixing  layer)  and  other  mechanisms  of  importance  in 
other  turbulent  flows  have  begun  to  be  recognized.  Although  it 
has  not  happened  yet,  these  insights  into  the  physics  of  turbu¬ 
lent  flow  are  expected  to  have  important  influences  on  turbulence 
modeling.  It  is  our  intent  to  contribute  to  the  understanding  of 
the  basic  mechanisms  underlying  turbulent  flows  by  making  use  of 
the  wealth  of  information  implicit  in  the  flow  fields  calculated 
by  full  numerical  simulation  of  the  Navier-Stokes  equations. 

Our  tools  in  this  study  are  primarily  visual:  by  making 
contour  plots  of  various  quantities  of  importance  in  these  flows 
(e.g.,  turbulent  stress  or  vorticity) ,  we  are  attempting  to  find 
evidence  of  a  mechanism  akin  to  pairing  to  explain  the  length- 
scale  growth  which  occurs  when  turbulence  is  strained  or  sheared. 
Specifically,  we  are  focusing  on  the  homogeneous  shear  flows 
calculated  by  Rogallo  (ref.  1) ,  treating  them  as  "experimental 
data"  measured  in  complete  detail  at  a  series  of  discrete  time 
steps. 
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In  the  past  year,  we  have  developed  the  data-handling  pro¬ 
cedures  and  computer-graphics  routines  necessary  to  investigate 
the  structure  of  homogeneous  shear  flows,  and  we  have  begun  this 
study  using  one  of  Rogallo's  calculated  cases  (ref.  1,  case  BSH9 , 
shear  rate  S  =  3u/9y  =  28) .  These  cases  were  computed  in  a 
coordinate  system  which  shears  with  the  mean  motion,  and  which 
is  remeshed  periodically  to  avoid  having  the  computational  volume 
become  too  small  in  any  direction  to  resolve  the  pertinent  turbu¬ 
lent  scales.  Velocity  fields  are  available  for  analysis  at  the 
instants  of  time  at  which  the  coordinate  system  is  orthogonal 
(St  =  0,1,2, •••).  Rogallo's  calculations  were  done  on  the  ILLIAC 

3 

computer  at  NASA/ARC  using  128  grid  points.  Our  data  reduction 
and  post-processing  is  done  on  the  ARC  CDC  7600;  restrictions  on 
storage  space  and  computer  costs  require  that  we  deal  with  a  64 
approximation  to  these  results,  which  is  obtained  by  retaining 
only  the  lower  half  of  the  wavenumbers  in  each  direction  resolved 
in  the  128 3  computation.  Obviously,  we  lose  details  of  the  fine 
structure  in  this  way,  but  the  large  scales  which  carry  the  bulk 
of  the  energy  and  which  are  responsible  for  most  of  turbulent 

3 

transport  are  retained  (90%  of  the  energy  in  the  128  field  is 

3 

retained  in  our  64  representation) . 

We  currently  can  produce  contour  plots  in  x-y,  x-z,  or  y-z 

2  2 

planes  of  each  component  of  the  turbulent  stress  (i.e.,  u  ,  v  , 

2  .  . 
w  ,  -uv,  -uw,  -vw)  ,  each  component  of  the  turbulent  vorticity  to., 

.  .  2  1 

and  the  square  of  the  magnitude  of  the  turbulent  vorticity  to  . 

In  our  investigation  into  the  structure  of  this  homogeneous  shear 
flow,  we  have  produced  several  planes  of  each  of  these  quantities 
at  several  instants  of  time  in  its  evolution. 

We  have  only  begun  to  make  use  of  these  "flow  visualization" 
techniques,  and  any  conclusions  we  could  draw  at  this  point  are 
very  preliminary.  However,  some  interesting  observations  can  be 
made,  which  are  illustrated  in  figures  8-12.  The  results  for  the 
shear  runs  in  these  figures  are  for  St  =  10  (near  the  end  of  the 
portion  of  the  simulation  of  this  case  which  Rogallo  claims  to 
be  adequately  resolved)  and  are  identified  with  the  label  SS1028 
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(St  =  10,  3u/3y  =  28).  In  these  plots,  contours  of  positive 
values  of  the  quantity  in  question  are  shown  with  solid  lines, 
negative  values  with  dashed  lines. 

Figures  8(a)  and  (b)  show  contours  in  an  x-y  and  x-z  plane, 
respectively,  of  the  predominant  turbulent  stress  in  this  flow, 

-uv.  The  diagonal  lines  of  figure  8(a)  are  drawn  parallel  to 
one  of  the  principal  axes  for  the  turbulent  stress  in  this  flow 
field  (calculated  to  lie  at  a  22°  angle  with  respect  to  the  x 
axis,  ref.  1).  These  figures  show  some  alignment  of  the  con¬ 
tours  with  this  principal  axis,  and  none  with  the  principal  axis 
of  the  mean  strain  which  is  at  45°.  Because  the  structures  are, 
in  the  mean,  inclined  in  x-y  planes  [fig.  8(a)],  they  appear 
somewhat  shorter  in  the  x-direction  in  x-z  planes  [fig.  8(b)]. 
Additionally,  the  structures  appear  somewhat  thinner  in  the  z- 
direction  than  in  the  direction  normal  to  the  principal  axis  in 
figure  8(a),  suggesting  that  they  resemble  irregular  "disks", 
with  their  thin  dimension  more  or  less  aligned  with  the  z  axis. 

The  relative  "plumpness"  of  the  structures  in  the  x-y  plane  is 
counter  to  intuition;  we  had  expected  the  structures  to  be 
sheared.  By  way  of  contrast,  contours  of  -uv  are  shown  in  an 
x-y  and  x-z  plane  for  an  isotropic  flow  in  figures  9(a)  and  (b) . 

In  this  isotropic  flow,  no  particular  structure  is  apparent. 

(Note  that  the  relative  number  of  contours  is  not  significant  in 
the  comparison  of  the  plots  for  these  two  flows  since  it  reflects 
the  arbitrary  choice  of  contour  levels  used.) 

Contours  of  -vw  in  an  x-y  plane  of  the  shear  flow  [fig.  10(a)] 
again  show  alignment  of  the  structures  with  the  principal  axis  for 
stress,  with  a  "thinness"  normal  to  this  axis  perhaps  more  in  keep¬ 
ing  with  one's  expectations.  The  x-z  plane  for  this  quantity 
[fig.  10(b)]  shows  no  particular  structure.  It  is  interesting  to 
note  that  the  average  value  of  -vw  in  this  flow  field  is  about  2% 
of  the  average  value  of  -uv.  Because  the  contours  of  figure  10 
are  at  the  same  intervals  as  those  of  figure  8,  one  can  infer  that 
the  average  lack  of  -vw  in  the  flow  field  is  not  due  to  absence 
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of  -vw  in  the  field,  but  rather  to  cancellation  of  positive  -vw 
by  an  almost  equal  abount  of  negative  -vw. 

Further  samples  of  the  results  generated  to  date  are  shown 
in  figures  11  and  12,  where  x-y  and  x-z  planes  of  the  x  component 

fig.  11)  and  the  z  component  fig.  12)  of  turbulent 

vorticity  are  shown.  Long  thin  structures  aligned  with  the 
principal  axis  in  the  x-y  plane  are  again  prevalent,  as  is  a 
lack  of  any  organization  in  the  x-z  plane.  Something  lacking  in 
these  plots  and  in  others  we  have  generated  but  have  not  included 
here  is  any  indication  of  vortex  pairing. 

At  this  point  in  our  study,  the  mechanisms  responsible  for 
the  production  of  turbulent  energy  and  the  growth  of  the  turbu¬ 
lence  length  scales  are  not  clear;  the  most  that  can  be  ventured 
is  that  these  mechanisms  seem  to  be  different  than  those  acting 
in  the  inhomogeneous  version  of  this  flow,  the  mixing  layer. 

The  elucidation  of  the  essential  similarities  and  differences 
in  these  flows  is  scheduled  to  be  continued  under  ONR  sponsorship 
in  another  program. 

4.  DEVELOPING  A  MODEL  OF  TURBULENCE  NEAR  A  WALL  FROM 
SOLUTIONS  OF  THE  NAVIER-STOKES  EQUATIONS 

Most  of  the  work  on  the  development  of  a  model  for  turbu¬ 
lence  near  a  wall  has  been  reported  in  an  interim  technical 
report  (ref.  6) .  The  work  will  be  briefly  summarized  herein. 

The  model  for  the  first  time  makes  direct  use  of  what  is  known 
of  the  structure  of  the  flow  in  the  vicinity  of  the  wall.  The 
model,  which  can  lead  to  improvements  in  both  Reynolds  averaged 
and  subgrid-scale  turbulence  models,  is  based  on  the  solution  of 
the  Navier-Stokes  equations  in  the  viscous  sublayer. 

The  major  goal  of  the  work  is  to  determine  velocity  boundary 
conditions  to  be  applied  at  the  edge  of  the  viscous  sublayer.  In 
the  work  reported  in  reference  6,  time  and  space  phase  relation¬ 
ships  in  the  edge  boundary  conditions  were  explored.  Only  first 
harmonics  were  included  in  the  work  reported  in  reference  6. 
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Subsequently,  preliminary  work  has  been  done  to  include  higher 
harmonics  of  the  velocity  fluctuations  at  the  edge  of  the  sub¬ 
layer.  The  first-harmonic  work  was  surprisingly  successful,  and 
shows  that  there  is  considerable  potential  for  coherent-structure 
modeling  as  a  means  of  computing  the  dynamics  of  viscous-sublayer 
turbulence . 


During  the  course  of  exploring  various  velocity  boundary 
conditions,  it  was  realized  that  two  separate  components  would 
be  required  to  simulate  realistically  the  dynamics  of  viscous 
sublayer  turbulence.  One  component  represents  the  small-scale 
eddies  (SSE)  responsible  for  the  principal  production  of  turbu¬ 
lence  and  Reynolds  stress.  The  second  component  represents 
organized  large-scale  eddies  (LSE)  which  interact  with  the  small- 
scale  eddies.  The  crux  of  this  coherent  eddy  model  is  the  con¬ 
struction  of  appropriate  time-  and  space-dependent  boundary 
conditions  for  the  three  fluctuating  velocity  components  ug+ , 
ve+,  we+  at  the  outer  edge  of  the  viscous  sublayer  (at  y  =  yg)  . 

A  simple  construction  has  been  used  to  date  for  purposes  of 
illustrating  the  potential  and  the  main  characteristics  of  this 
type  of  turbulence  modeling.  With  all  quantities  expressed  in 
conventional  dimensionless  wall  variables,  the  velocity  boundary 
conditions  developed  are 


Component  1  —  SSE  Component  2  —  LSE 


u 


e+ 


=  2a, 


sin 


NjT 


sin  5 


+  /2(a^  -  a?)sin(N 


u2 


T  + 


*u2> 
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e+ 


we+ 


-28  sin  N^T  sin  c 
28  cos  N^T  cos  c,  +  Jl  (y2 


6^) sin 


N 


u2 


T  +  <j> 


w2 


(12) 


where  a,  8,  and  y  are  the  rms  intensities  of  fluctuation  for  ug+, 
vg+,  we+/  respectively,  is  the  mean  frequency  of  SSE  burst 
events  (ejection/sweep  events),  N  £  is  the  mean  frequency  of  LSE, 
C  =  2ttZ/X  is  the  dimensionless  spanwise  distance  variable,  T  is 
the  dimensionless  time,  X  is  the  mean  spacing  between  high-speed 
and  low-speed  streaks,  a^/a  is  the  uv  correlation  coefficient 
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(0.45)  at  yg/  and  4>u2  and  ^w2  are  phase  angles.  The  flow  model 
is  periodic  both  in  time  and  span.  The  values  of  the  various 
constants  used  in  the  comparisons  presented  in  reference  6  are 
listed  in  table  4. 


The  computational  mesh  used  for  calculations  contained  17 
equally-spaced  points  in  the  spanwise  direction  and  24  points  in 
the  vertical  direction  clustered  near  the  wall  and  near  the  edge 
by  the  functions 


and 


where 


for 


o  1  . 

yi  =  2< 


r 


1  +  —  tanh  £ .  tanh 
a  |_;i 


o 

,  =  !iy 
1  e 


£  =  2  (i-  ~1  -P.  _  i 

H  N  +  3 


(13) 


i  —  1/2,3, 


/N 


and 


a  =  0.98 

This  construction  of  velocity  boundary  conditions  was  guided 
by  experimental  observations  of  organized  sublayer  structure  made 
during  the  past  two  decades.  For  example,  the  sin  z,  factor  for 
ufi+  corresponds  to  the  observation  of  high-  and  low-speed  streaks 
spaced  spanwise  a  mean  distance  of  X  apart,  while  the  factors 
sin  c,  -sin  £ ,  and  cos  £  in  the  ug+,  ve+/  and  wg+  equations, 
respectively,  correspond  to  a  simple  modeling  of  the  observations 
of  contra-rotating  vortical  motion  near  a  wall.  The  absence  of 
a  phase  angle  in  time  between  Ue+  and  vg+  corresponds  to  the 
observation  that  the  uv  correlation  in  the  log  region,  and  hence 
at  y  =  yg,  is  maximum  with  zero  time  delay  between  their  respec¬ 
tive  signals.  The  180°  phase  difference  between  ug+  and  vg+ 
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corresponds  to  observations  from  conditional  samples  in  the 
sublayer  that  u  and  v  are  180°  out  of  phase  during  the  Reynolds- 
stress  intensive  ejection-sweep  event.  The  90°  phase  difference 
in  time  between  vg  and  wg  corresponds  to  the  observation  that 

~2 

the  derivative  Ovg/3y)  at  the  outer  edge  of  the  viscous  sub¬ 
layer  is  zero.  A  comparison  of  a  number  of  turbulence  charac¬ 
teristics  computed  from  this  simple  model  with  those  measured  in 
experiments  showed  surprisingly  good  agreement  at  the  Reynolds 
number  investigated  (Re^  ~  29,000). 

It  is  of  interest  to  note  that  the  two-component  type  of 
velocity  boundary  condition  at  the  outer  edge  of  the  viscous 
sublayer  is  precisely  compatible  with  the  concept  of  "active" 
and  "inactive"  components  of  turbulent  motion  in  the  log  region 
as  characterized  by  Townsend  (ref.  18)  and  Bradshaw  (ref.  19). 
Component  1  is  active,  involving  small-scale  eddies,  producing 
the  Reynolds  stress,  being  rotational,  and  being  dependent  on 
wall  variables.  Component  2  is  inactive,  involving  large-scale 
eddies,  producing  energy  but  no  Reynolds  stress,  being  irrota- 
tional,  and  being  dependent  on  outer  variables.  There  is  much 
experimental  information  that  turbulent  flow  comprises  these 
two  distinct  types  of  motion  in  the  log  region,  and  hence  at  the 
outer  edge  of  the  viscous  sublayer  where  the  boundary  conditions 
are  applied. 

Some  preliminary  runs  were  made  for  Re^  =  11,000,  as  well 
as  for  Re^  =  29,000.  The  results  were  not  significantly  differ¬ 
ent  for  this  relatively  small  range  of  Re^,  and  only  the  compu¬ 
tations  for  Re^  =  29,000  were  presented  in  reference  6.  These 
latter  computations  correspond  essentially  to  the  Re^  of  Lauf er ' s 
classical  pipe-flow  experiments  (ref.  20),  and  the  Re^  near  which 
many  other  viscous-sublayer  experiments  were  conducted.  It  is 
believed  that  computations  over  a  wide  range  of  Refi  should  be 
delayed  until  the  intermittent  production  of  Reynolds  stress  can 
be  properly  built  into  the  velocity  boundary  conditions.  Compu¬ 
tations  at  high  Re^  will  require  much  larger  amounts  of  computer 
time.  A  run  for  Re,  =  29,000  takes  about  4-5  minutes  on  the 


CDC  7600  computer,  while  a  run  for  Re^  =  290,000  will  take 
roughly  45  minutes.  Hence  it  is  prudent  to  essentially  complete 
the  model  development  before  runs  at  very  high  Re^  are  made. 

Computations  at  Re^  =  29,000  were  made  for  two  different 
pressure  gradients  corresponding  to  pipe  flow  and  channel  flow. 
The  results,  presented  in  reference  6,  showed  that,  for  the  mild 
pressure  gradients  of  these  two  flows,  the  dynamics  of  viscous 
sublayer  turbulence  are  nearly  indistinguishable  from  those  for 
zero  pressure  gradient  on  a  flat  plate.  These  results  are  con¬ 
sistent  with  experimental  observations  for  which  measuring 
inaccuracies  have  masked  the  very  small  differences  in  viscous 
sublayer  flow  which  exist  between  pipes,  channels,  and  flat 
plates . 

In  work  done  subsequent  to  the  interim  report  (ref.  6),  the 
computer  program  was  modified  to  compute  the  spatial  derivatives 
of  velocity  that  determine  the  rate  of  energy  dissipation.  No 
results  were  obtained  to  evaluate  the  ability  of  the  model  to 
calculate  the  dissipation. 

The  primary  weakness  of  the  model  reported  in  reference  6 , 
as  mentioned  briefly  earlier,  is  that  it  does  not  simulate  the 
known  physics  of  intermittent  Reynolds  stress  production.  The 
initial  model  corresponds  to  Reynolds  stress  being  produced  in 
a  simple  sinusoidal  fashion  with  time;  whereas  in  reality  the 
viscous  sublayer  is  quiescent  most  of  the  time  and  almost  all 
of  the  Reynolds  stress  is  produced  in  about  10-20  percent  of  the 
time  through  energetic  burst/sweep  events.  Thus  the  main  objec¬ 
tive  of  work  to  improve  the  model  is  to  develop  modified  velocity 
boundary  conditions  that  simulate  realistically  the  observed 
intermittent  nature  of  Reynolds  stress  production. 

The  concept  for  accomplishing  this  objective  is  relatively 
simple.  Reynolds  stress  is  produced  only  by  the  active  component 
of  small-scale  eddies  (SSE) ;  hence  only  the  time  function  sin  N^T 
in  the  boundary  condition  for  ug+  and  vg+  (and  perhaps  also  the 
cos  N^T  function  for  wg+)  will  need  to  be  modified.  The  proposed 
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method  of  modification  is  illustrated  in  figure  13.  Additional 
harmonics  of  the  Fourier  series  that  approximates  a  step-function 
burst/sweep  event  have  been  added  to  the  sin  N^T  term  in  the 
velocity  boundary  conditions  in  order  to  simulate  intermittent 
production  of  Reynolds  stress.  Studies  of  the  time  and  spatial 
|  phase  relationships  of  the  new  boundary  conditions  are  scheduled 

to  be  conducted  under  ONR  sponsorship  in  another  program. 

5.  CONCLUDING  REMARKS 

Studies  have  been  conducted  on  several  aspects  of  the 
structure  and  computation  of  incompressible  turbulence.  In  the 
first  of  these,  involving  the  testing  in  strained  and  sheared 
turbulence  of  subgrid-scale  models  for  use  in  the  large-eddy 
simulation  of  homogeneous  turbulent  flows,  it  has  been  shown 
that  eddy  viscosity  models  are  quite  inaccurate  when  severe 
strain  has  been  imposed  on  the  flow.  If  the  model  is  to  be  used 
in  strained  or  sheared  turbulence  despite  this,  it  is  better  to 
remove  the  mean  strain  or  shear  from  the  velocity  field  before 
the  model  is  computed.  The  deficiency  in  eddy  viscosity  models 
appears  to  come  from  the  fact  that  the  principal  axes  of  the 
exact  subgrid-scale  stress  and  the  modeled  values  are  not  aligned. 
New  models  based  on  the  notion  of  scale  similarity  have  been 
tested  and,  when  used  in  combination  with  other  models  to  provide 
energy  dissipation,  appear  to  offer  considerable  promise.  And 
finally,  it  has  been  observed  that  when  the  strain  is  strong,  the 
energy  cascade  is  from  the  small  scales  to  the  larger  ones,  a 
result  with  possible  considerable  importance  for  turbulence 
modeling . 

The  second  study  involves  the  search  for  the  mechanism  under¬ 
lying  the  growth  of  scales  in  homogeneous  shear  flow.  Data 
management  and  computer-graphic  routines  have  been  developed  to 
allow  "flow-visualization-like"  examination  of  flow  fields  cal¬ 
culated  by  full  numerical  simulation  of  the  Navier-Stokes 
equations.  Preliminary  work  accomplished  to  date  in  this  study 
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suggests  that  the  mechanisms  controlling  homogeneous  shear  flow 
are  different  than  those  underlying  the  inhomogeneous  version 
of  this  flow,  the  mixing  layer. 

The  third  study,  that  of  flow  in  the  viscous  sublayer 
revealed  a  good  potential  for  constructing  quantitative  models 
of  viscous  sublayer  turbulence  founded  on  experimental  observa¬ 
tions  of  the  physics  of  organized  quasi-periodic  eddy  structures. 
Much  physical  information  has  been  accumulated  in  the  past  two 
decades  on  organized  sublayer  structures.  A  sizeable  part,  but 
not  all  of  the  essential  physics,  has  been  included  in  the 
initial  model.  Embodiment  of  physics  of  coherent  structures, 
however,  has  not  been  possible  within  the  framework  of  conven¬ 
tional  Reynolds-average  modeling.  It  is  clear  that  the  present 
approach  to  modeling  can  be  applied  to  study,  for  example,  the 
effects  of  Reynolds  number  on  turbulence  intensities,  dissipation, 
etc.,  as  well  as  the  effects  of  pressure  gradients  much  stronger 
than  those  encountered  in  channel  and  pipe  flow.  Moreover,  the 
method  can  be  extended  to  more  general  conditions,  such  as  to 
flows  with  heat  transfer  and  compressibility.  Modeling  the  ef¬ 
fects  of  streamwise  surface  curvature,  however,  would  apparently 
require  extension  to  full  three-dimensional  flow  by  including 
the  streamwise  velocity  derivatives  neglected  in  the  initial 
model . 
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TABLE  1.-  ENERGY  TRANSFER  TO  THE  SUBGRID-SCALES 

BY  T . .  AND  X . .  STRAIN  RATIO  =  4 . 

13  ij  . 


Grid 

T  (sec-1) 

3u .  2  - 

<Tij  9x  j >  (CIT1  sec  } 

3u  •  2 

<x.  .  (cn\  sec 

ID  3x^ 

643 

5.30 

-2.85 

3.62 

643 

10.61 

-4.77 

-0.18 

1283 

10.61 

-8.30 

-4.63 

TABLE  2.-  CORRELATION  COEFFICIENTS,  ISOTROPIC 

DECAY  (R  =143). 
sgs 

Smagorinsky  Equation  (10) 


Tensor 

0.30 

0.52 

Vector 

0.27 

0.46 

Scalar 

0.37 

0.45 

TABLE  3.-  CORRELATION  COEFFICIENTS,  PLANE  STRAIN 

(ert  =  4.01,  Tt,  =  3.75) . 
a 


Tensor 


Smagorinsky  Equation  (10) 


0.11 


0.49 


Vector  0.11 


0.47 


Scalar 


-0.09 


0.45 


TABLE  4.-  VALUES  OF  CONSTANTS  USED  IN  INITIAL  MODEL  OF 
TURBULENCE  IN  THE  VISCOUS  SUBLAYER. 

a  =  2.0 

=  0.9 

6  =  1.0 

Y  =  1.3 

Nx  =  .025 

N  ,  =  .025 
u2 

*u2  =  77/3 
*w2  =  371/4 


kl 
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STRAIN  RATIO 


Figure  1.-  Energy  transfer  via  backscatter 
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Figure  2.-  Scalar  correlations  for  Sirtagorinsky  model. 
Open  symbols:  mean  strain  rate  included  in  model. 
Closed  symbols:  no  mean  strain  rate  in  model. 
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Figure  3.-  Scalar-level  parameter,  Smagorinsky 
model.  Mean  strain  rate  included  in  model. 
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Figure  4.-  Scalar-level  parameter,  Smgorinsky 
model.  No  mean  strain  rate  in  model. 
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Figure  8(b) 
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